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Abstract 

Let A: be a differential field with algebraic closure k, and let [A] : 
Y' = AY with A 6 A4 n (k) be a linear differential system. Denote by 
the Lie algebra of the differential Galois group of [A]. We say that 
a matrix R G A4 n (k) is a reduced form of [A] if R G g(fc) and there 
exists P G GL n (k) such that R = P _1 (AP - P') G fl(fc). Such a form 
is often the sparsest possible attainable through gauge transformations 
without introducing new transcendents. In this article, we discuss how 
to compute reduced forms of some symplectic differential systems, arising 
as variational equations of hamiltonian systems. We use this to give an 
effective form of the Morales-Ramis theorem on (non-)-integrability of 
Hamiltonian systems. 



1 Introduction 

This article lies at the crossroads of differential Galois theory and the complete 
(Liouville) integrability of Hamiltonian systems. Let [A] : Y' = AY with 
A G A4 n (k) be a linear differential system, where k denotes a differential field 
of characteristic zero with algebraically closed constant field C . 

On the differential Galois theory side, we propose (following works of Kolchin, 
Kovacic, Singer, Mitschi and others) a notion of reduced form for the system 
[A]. Let G denote the differential Galois group of Y' = AY and g its Lie al- 
gebra. We say that a matrix R G A4 n (k) is a reduced form of [A] if R G g(k) 
and there exists P G GL n (k) such that R = P^ 1 (AP — P'). Such a form turns 
out to be very natural and somehow the most concise attainable through gauge 
transformations without introducing new transcendents. When a system is in 
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Hartmannl (120051); I Juan and Ledetl (|200/t . see more 



van der Put and Singerl (|2003j) , generally for inverse 



reduced form, many of its intrinsic properties can be readily computed. In our 
case, we are mainly concerned with detecting the (non-)abelianity of the Lie al- 
gebra of the differential Gal ois group of \A]. Stemm ing from works by Kovacic 
and Kolchin ( Kolchi nl (|l999h or lKovacid on the inverse problem of 

differential Galois theo ry, th is notion of reduced form has been since developed 
bv lMitschi and Singerl (Il 996l. l2002h: ICook et al.1 (120051) ffor inv erse p roblems as 
well) and m any oth e rs (e. g 
references in lSinger (l2009h : 
problems). We explore here both some of the advantages and the constructibil- 
ity of this reduced form applied to the context of direct problems in differential 
Galois theory. 

On the Hamiltonian system side, the Morales-Ramis theorem states that if a 
Hamiltonian system is meromorphically completely (Liouville) integrable, then 
the Lie algebra of its variational equation (a linear differential system) along a 
non-constant integral curve is abelian. Applying the above notion of reduced 
form (effectively), we are able to decide whether this Lie algebra is abelian or 
not, and thus propose new (non-)integrability criteria. When the Lie algebra is 
abelian, our reduced form is useful to simplify the study of higher variational 
eq uations in order to apply the Mora l es-Ramis-Simo theo r em (t his is explored 
in lAparicio- Monforte and Weill (|201ll) ; lAparicio-Monfortej (|2010l) ). 

Many of the existing examples are Hamiltonian systems with two degrees of 
freedom. In such case, we will show (or recall to specialists) how the normal 
variational equation can be put into a reduced form using Kovacic's algorithm. 
In the literature, many authors take into consideration only the normal varia- 
tional equation [N]. Since the Lie algebra of [N] is only a quotient of g, the 
data thus obtained cannot be complete, especially as far as the abelianity of g 
is concerned. 

In this paper, we give an algorithm which, by looking for a reduced form of a 
linear differential system, will either show that the full system has non-abelian 
Lie algebra (hence proving non-intcgrability of the original differential system) 
or (in the abelian case) return a reduced form for the variational equation. 

Our construction is quite systematic an d should generalize to many o t her sy s- 
tems (such generalizations are initiated in lAparicio-Monforte and Weill (|201lf) ). 
It also shows the lovely simple structure of reduced forms. This will be very 
useful in studying higher variational equations. Systems occurring in higher vari- 
ational equations are reducible and of big sizes; computi ng their Galois groups 
is difficult. Although a general decisi on procedur e exist s (|Compoint and Singer 
(1999) in completely reducible cases. iHrushovski ( 2002 ) in general), studies to- 
wards extensive description s of dif ferential Galois groups in reducible cases (e.g 
Bertrandl (l200lh: iHardouinl (l2005l)1 are still in pro g ress. The case of two com- 



Bermanl (|2002f )) is the only 



pletely reducible factors (jBerman and Singer (1999); 
one which seems to be fully understood. We propose an approach which focuses 
on finding abelianity criteria for the Lie algebra instead of computing the Galois 
group. 

This paper is split into four sections apart from this introduction. Section 
2 deals with theoretical background material such as differential Galois theory 
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and integrability of Hamiltonian systems. Section 3 introduces the concept 
of reduced form as well as its importance and usefulness, showing that the 
Kovacic algorithm can be re-used to put second order systems into a reduced 
form. In Section 4, we consider linear differential systems whose matrix A lies 
in sp(4, k) (typically the variational equation of a Hamiltonian system with 
two degrees of freedom) and provide a reduction algorithm together with an 
abelianity criterion. In section 5, we show how to apply our techniques to 
questions of integrability of Hamiltonian systems. In particular, our method 
allows us to reprove that the lunar Hill system is not meromorphically integrable. 
Acknowledgements: we warmly thank S. Simon for fruitful conversations and 
for suggesting the study of the Hill example which is given in section 5. We also 
thank Thomas Cluzeau, Elie Compoint, Maria Przybylska, Andrzej Maciejewski 
and particularly the referees for their useful suggestions. 

2 Background material 

This section contains necessary background material and no new results. 



2.1 Some differential Galois theory 

General references for this section are Ivan der Put and Singer (2003); Singer 



(2009) and many others, cited therein. Let ( k , ' ) be a differential field with 



an algebraically closed field of constants C of characteristic zero. Let [A] 
Y' = AY denote a linear differential system with A £ M n (k). A Picard-Vessiot 
field K is a minimal differential extension of k generated by the entries of a 
fundamental solution matrix U of [A] . The differential Galois group of A is the 
group of the automorphisms over K that leave k invariant and that commute 
with the derivation: 

The Galois group G is a linear algebraic group acting on the vector space of 
solutions of [A] and as such admits a faithful representation in GL n (C) once a 
fundamental solution matrix U is chosen. The Lie algebra of G is q := Lie(G) = 
T e G, a C-vector space of matrices endowed with a Lie algebra structure by the 
usual Lie bracket [A, B] := AB — BA. The Lie algebra g is abelian if and only 
if the connected component of the identity G° is abelian as well. 

We say that two systems Y' = AY and Z' = BZ with A, B £ M n (k) are 
gauge equivalent if there is a linear change of variable Y = PZ with P £ GL n {k) 
such that Z' = BZ, i.e: 

B = P[A] := P~ 1 (AP - P') (gauge transformation). 

Two gauge equivalent systems share the same Galois group G. Moreover, V :— 
P^ 1 U is a fundamental solution matrix of [B] ; given a £ G with matrix M a , 
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we have a(U) = UM a and a^ 1 ) = P" 1 so a(V) = P- x UM a = VM a , i.e 
the representation of G is unchanged on this new system. If the coefficients of 
the gauge transformation belong to the algebraic closure k, then G is altered 
whereas g and G° remain unchanged. 



2.2 Hamiltonian Systems 

Let (M, w) be a complex analytic symplectic manifold of complex dimension 
2n. By the Darboux theorem, we know that, locally, M is isomorphic to an 
open domain U C C 2n and that, taking locally suitable coordinates (q,p) = 
(qi, . . . , q n ,Pi, ■ ■ ■ ,Pn), we can deal with associated entities (such as Hamiltons 
equations and the Poisson bracket) in terms of the matrix 



J 



On In 

In. 



In these coordinates, given a function H G C 2 (U) : U C, we define a 
Hamiltonian system over U C C 2n as the differential system given by the vector 
field X H = J ■ VH: 

qi = ^~{q,p) , Pi = ~'jr-{q,p) for i = l...n. (1) 

Opt dqi 

Let z(t) be a parametrization of an integral curve T C U that satisfies ((T|). The 
Hamiltonian H is constant over those integral curves. Indeed, 

X H ■ H := (Vi? , X H ) = (VH , JVif) = 0. 

Therefore integral curves will lie on the levels of energy of H . A function 
F : U — > C, meromorphic over U, is called a meromorphic first integral of 
(p} if Xh ■ F = (i.e. F is constant over the integral curves of H). The 
Poisson bracket { , } of two meromorphic functions /, g defined over a symplectic 
manifold is defined as {/,<?} '■— (Xf,S7g) = (—X g ,Vf) and, in coordinates, 



% <9pi <9pi dqi 



The Poisson bracket endows the set of first integrals of ([T]) with a Lie algebra 
structure ; in fact, a function F is a first integral of ([T]) if and only if it is in 
involution with H , i.e. {F, H} = 0. 

A Hamiltonian system is called meromorphically Liouville integrable if it 
possesses n first integrals Hi = H , . . . , H n meromorphic over U satisfying: 

- they are functionally independent: V-ffj., . . . , VUn are linearly indepen- 
dent over U, 

- they are in involution: {Hi , Hj} = for i,j = l...n. 
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The theory of Mo r ales and Ramis, d evelo pi ng on founding works of Ziglin. 
and followers (IZiglinl (Il982l Il983lh jltol (Il985lj. I Yoshidal (119861 Il987bl lal. Il988h 



Baider et all |l996); Ch urchill et al.l (|l995h . iMorales Ruizl(j 19991) ) . aims at prov 



ing rigorously non-integrability using differential Galois groups (or monodromy 
groups) of variational equations. 



2.2.1 Variational equations 

Let r C U be an integral curve of ([TJ parametrized by z(t). The differential 
field k := <C(z(t)} will be called the coefficient field (or, informally, the field 
of rational functions). We define the variational equation of ([TJ) along T as 
the linearization of (JTJ) along z(t). ft describes the behavior of the solutions 
of fTJ) near z(t). In other words, if zo(t) and z\ (t) := zo(t) + Y(t) where Y{i) 
is an infinitesimal perturbation of zo(t), then a first order approximation gives 
Y' = d Zo ( t )XHY(t). Taking coordinates, the variational equation along z(t) is 
given by 

Y' = AY with A := J ■ Ress(H)(z(t)). (3) 

As the Hessian Hess(H)(z(t)) is a symmetric matrix, we have A € sp(2n, k) 
(i.e A T ■ J + J ■ A — 0). Thus, as recalled in the Appendix, there exists a 
fundamental matrix of solutions U of Y' = AY in a Picard-Vessiot field K 
such that U £ Sp(2rt, K). In the sequel, we will work with such a symplcctic 
fundamental solution matrix U(t) G Sp(2n,if). 

We denote by G the differential Galois group of the variational equation ^ 
and by g its Lie algebra. As system ([3]) is Hamiltonian (A £ sp(2ro, fc)), G is a 
subgroup of Sp(2n, C) and g C sp(2n, C) (see next section). 



Theorem 1 (Morales and Ramis, see Morales Ruid (|l999l ) Theorem 4.2 p.81 



Let z(t) be a non-singular integral curve of the Hamiltonian system {!]). Let 
(0) be its variational equation along z(t). If (Qp is meromorphically Liouville 
integrable then the Lie algebra g of (pj) is abelian. 

This theorem is in fact a non-integrability criterion: when checking whether fj 
is abelian only negative answers are conclusive. Indeed, there are non-integrable 
systems for which the Lie algebra (of the differential Galois group) of their varia- 
tional equation along a given integral curve is abelian In this case, one is bound 



to consider hi gher order variational equations, see e.g Morales- Ruiz et al. (2007) 
section 3.4 or lAparicio-Monforte and Weill (|201lh : lAparicio-Monfortel (|2010h 



2.2.2 The Normal variational equation 

In general, studying the abelianity of g is simpler than finding either a complete 
set of meromorphic first integrals of ([TJ or an obstruction to their existence 
without taking into account (j3]). Still, checking straightforwardly the abelianity 
of g is not an easy matter. 

However, we can take advantage of the fact that z'(t) is a particular solu- 
tion of the variational equation ([3]) along z(t). Indeed, Proposition 4.2 p 76 in 
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Morales Ruiz ( 19991) (see also our Appendix) ensures the existence of a sym- 
plectic gauge transformation that allows us to reduce this variational equation 
(i.e. to rule out one degree of freedom) and to obtain the normal variational 
equation (NVE). In the new coordinates, (NVE) can be written as Z' = NZ 
where N G sp(2(n — 1); k): (NVE) is therefore yet another (Hamiltonian) linear 
differential system. 

Consider now Un S Sp(2(n — 1),Kn) a fundamental matrix of solutions of 
Z' = NZ, where Km D k is a Picard Vessiot extension for (NVE). We have 
k C K N C K so G := dAut KN {K) <\ G := dAut k (K). Since the differential 
Galois group of the (NVE) is given by Gm '■= dAutk(Kjy) ~ G/G then the Lie 
Algebra of Gn is Qn = fl/fl- Thus, if qn is non-abelian then g is not abelian 
either. In fact, the usual method, ever since Morales introduced it, consists in 
reducing the ([3]) to its (NVE) and then argue about the abelianity oi Qn- If Qn 
turns out to be abelian, no conclusion for §5§ is reached and so the higher order 
variational equations are explored. We will improve this procedure in sections 
3 and 4. 

Let us illustrate the notion of normal variational equation by taking a general 
example for 2n — 4. In such a case, the variational equation ([3]) can be written: 



Y' = AY with A := 



an 
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esp(A,k). 



(4) 



Take a particular solution z — (z\, Z2, 23, z^) T of ([1} (for n = 2). Then 



(A 



1 2 ' 3 ' 



z' 4 ) T is a particular solution of ((4]). Completing z' to a symplcctic 



basis of k 4 (see Appendix) yields the symplectic change of variable Z = PU, 
where 



P 
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eS P (A,k). 



(5) 



It induces a (symplectic) gauge transformation A^ 
sp(4; fc), where 



:= P[A] = P- 1 {AP-P') e 



A N := 






ai2,n 


Cl3,n 
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(6) 



The algebraic normal variational (or for convenience normal variational 
equation, see Remark equation is given by t/ = where 



N 



n 2 i 



ni2 



G 



3 Reduced form of a linear differential system 



3.1 A Kovacic reduced form 



The study of the integrability of dynamical systems has origi nated numerous 
deep studies on local normal forms near equilibrium points (e.g Birkhoff (1966) 
and references therein) or even along periodic solutions (e. g. iKozlovl (|l996l) ). 
However, we do not know of a similar global notion properly defined along 
non-constant regular solutions. 

In this work, we propose to explore a weaker global notion, a notion of 



reduced form of a linear diffe rential system which is strong l y insp ired by the 
work of Kovacic a nd Kolchin ( Kovacic (Il969l 1971D , lKolchinl (|l999h l and more 
recent works like iMitschi and Singerl ( 19961120021) on the inverse problem in 



differential Galois theory. 

Definition 2. Let A belong to M n (k). Let G be the differential Galois group of 
Y = AY and g its Lie algebra. We say that A is in reduced form if A G g(k). 

The expression A g g(k) is read "A is a fc-point of $j" . The Lie algebra g is a 
vector space of dimension d spanned by a set of matrices Mi, . . . , M d C .M„(C). 
The set g(k) of k-points of g is defined as 

g{k) := {fxM x + ■■■ + f d M d , f t G k}. 

Similarly, given a linear algebraic group G, the set G{k) of fc-points of G is 
the set of matrices whose entries are in k and satisfy the defining equations of 
G. The question is whether a reduced form exists. If k is a Ci-fieldJ then the 
following result due to Kolchin and Kovacic shows that the answer is positive: 

Theorem 3 (see lKovacicl (|l97ll ) or Ivan der Put and Singerl (|2003l ) p.25 Corol- 
lary 1.32). Let k be a differential C\-field. Let A G Ai n (k) and assume that the 
differential Galois Group G of the system Y' = AY is connected. Let g be the 
Lie algebra of G. Let H be a connected algebraic group such that its Lie alge- 
bra f) satisfies A G \)(k). Then G C H and there exists P G H{k) such that the 
equivalent differential equation F' = RF , with Y — PF and R = P^ 1 (AP — P'), 
satisfies R G g(k). 

Corollary 4. We use the same notations as above but now assume that G is 
not connected. Let K be a Picard-Vessiot extension of k for Y' = AY. Let k\ 
denote the algebraic closure ofk in K . Then there exists P G H{k\) such that the 
equivalent differential equation f — Rf with Y = Pf and R := P^ 1 (AP — P') 
satisfies R G g(fci). 

Proof. By virtue of van der Put and Singerl ( 20031 ) (Proposition 1.34 p. 26) the 
connected component of G containing the identity G° , satisfies G° — dAutk ± (K). 



A field k is called quasi-algebraically closed (or Ci) if every non-constant homogeneo us 
polynomial P G k[X\ , . . . , X n ] of degree less than n has a non-trivial zero in fc n llLanel jl952f) V 
An instance of one such field is k = C(x). In addition, by Theorem 5 of Lang (1952) wc know 
that the algebraic closure of a Ci field is Ci as well. 
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As ki is an algebraic extension of k, it is still a Ci-field. We now pick k\ as 
a base field; then K is a Picard Vessiot extension of k\ with Galois group G° 
satisfying the hypotheses of Theorem [3] □ 

Thus, for a Ci-field k and A G A4 n (k), a reduced form of [A] will be given by 
a gauge equivalent matrix R £ g(fci), i.e R := P[A] for some P £ GL n (k\). The 
change of variable given by P is algebraic ; the related gauge transformation 
does not introduce transcendental coefficients, hence preserving G° and g. 

Definition 5. Consider a differential system Y' — AY . We keep the notations 
of Theorem^ and corollary^ A matrix P £ GL n (k,i) is called a reduction 
matrix for [A] if P[A] is in reduced form, i.e P[A] € 0(fci). 

In this work, we will show how to make the Kolchin-Kovacic reduction the- 
orem effective for 2x2 systems with Galois group in SL(2,C) and for some 
systems with group in Sp(4, C) (obtained via the linearization of Hamiltonian 
systems). In both cases, we will see that one does not need the "Ci -field" hy- 
pothesis. There are other types of systems (e.g 3x3 systems with Galois group 
SO (3)) for which effective reduction may be per formed but where the "Ci -field" 
hypothesis cannot be eluded (see chapter 3 of Aparicio-Monfortel (2010), also 
lJuan and Ledetl l|2007l )). 



3.2 The Lie algebra associated to A 

Consider a matrix A = G M n (k). Let ai, . . . ,a r denote a basis of the 

C-vector space spanned by the a^j for i,j = 1, . . . ,n. We thus have a decom- 
position A :— J2i=i a i(t)Mi where the Mj are constant matrices. Of course, 
the choice of the dj in this decomposition is not unique, but the C-vector space 
generated by the Mi is. 

The Lie algebra generated by Mi, . . . , M r (i.e generated as a C-vector space by 
the Mi and all their iterated Lie brackets) will be called the Lie algebra associ- 
ated to A(t) and denoted by Lie(A). Its dimension d (as a vector space) satisfies 
r < d < n 2 . With this terminology, a system Y' = AY is in reduced form if 
Lie(A) is the Lie algebra Lie(Y' = AY) of the differential Galois group. 
This notion of Lie algebra associated to A appe ars in works of Magnus and of 
Wei and Norman, IWei and Normanl (|l963l . Il964 . 



Pick a gauge transformation P e GL„(fc) and let P[A] := P~ l {AP — P'); we 
say that P[A] is a partial reduction (and that P is a partial reduction matrix) 
if Lie(P[A}) C Lie(A). As we will see in section 4, a reduced form is obtained 
by means of a sequence of partial reductions. 

Consider now a system Y' = RY in reduced form, i.e R £ g(k) where 
g := Lic(F' = RY). Call r a minimal number of generators of g as a Lie 
algebra, (which is generally less than its dimension d as a vector space) and 
let m be the dimension of the C-vector space generated by the coefficients of 
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R. We say R is maximally reduced if m = r. Note that if q is abelian, then a 
reduced form is automatically maximally reduced (because r = d). 

Example 6. Consider the linear differential system Y' — AY with 



A := 



fi 





f3 
f-2 



-fl 



e Mi{k) 



(where functions fi S k are linearly independent over C) and assume that this 
is a reduced form, i.e the Lie algebra of its Galois group is the Lie algebra 
generated by 



Mi 



M 2 



and M 3 



We see that [Mi,M 2 ] = 2M 3 , which is linearly independent from Mi and M 2 . 
Thus, if we call g the Lie algebra generated by Mi and M2, it is spanned over 
C by Mi, M 2 and M 3 and we have that dimc(g) = 3 > 2. We see that system 
[A] is maximally reduced if (and only if) / 3 = O.D 



3.3 Solving abelian reduced systems 

Consider a linear differential system Y' = M(t)Y such that M(t) € f](fc) where 
f) is the Lie algebra associated to M(t). If \) := span c {Mi, . . . , Mj} then we 
can write 



M 



i=i 



fiMi with fi e k for i = 1 . . . d. 



As noted by Wei and Norman (who use this in Wei and Norman ( 1963 , 19641 ) ) , 
if t) is abelian, solving the system is straightforward. For instance, suppose 
that d = 2 and that f) is abelian, i.e. [Mi , M2] = 0. In this situation the 
system will be of the form Y' = {f\M\ + faM-zjY ' . The two differential sys- 
tems: Y{ = fiM{Yx and F 2 ' = f 2 M 2 Y 2 admit respectively Ui := exp(//iMi) 
and U2 := exp(J f 2 M 2 ) as fundamental solution matrices. As [U\ , M 2 ] = 0, 
a calculation shows that {U\U 2 )' = (/1M1 + f 2 M 2 )UiU 2 . This argument can 
be extended to d > 2 by induction. Therefore, if t) is abelian and we con- 
sider a system M in reduced form, we only need to solve each separate system 
Yj = fjMjYj-. U — nf=i ex P(J fiMi) is a fundamental solution matrix for the 
complete system. Such solving methods apply, more generally, when A com- 
mutes wit h f A or when Li e (A) is solva ble (this is the Wei-Norman method 
exposed in IWei and Normanl (|l963lll964) V 

Of course, solving is not the ultimate aim of system reduction, but these 
formulae will be useful in proving the correctness of our reduction procedure. 



3.4 Reduced form when G C SL(2, C) 



The Kovacic algorithm iKovacid (|1986l ) is an algorithm devised to compute Li- 
ouvillian solutions of second order linear differential equations. In fact, it can 
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be used to compute reduced forms for Y' = NY with N E A4 2 (k), when 
G° N C SL(2,C). We recall standard notations for abelian (non-trivial) con- 
nected subgroups of SL(2,C): the additive group is G a := { [ J \ ] : c G C} 
with Lie algebra Q a := { [ £ J ] : aeC}, and the multiplicative group is 

G m := { [ o c -i ] : c G C*} with Lie algebra m := {[ 2 _° Q ] : a G C}. 

In what follows, we assume that the differential Galois group of [N] is in 
SL(2,C) and Tr(iV) = 0, i.e N E sl(2,k) (otherwise an easy reduction puts 
us in this form). As a consequence, any fundamental solution matrix has a 
constant determinant. 

Let Km be a Picard-Vessiot extension of k associated to Y 1 = NY. The next 
proposition gives a complete reduction procedure, using the Kovacic algorithm, 
i.e shows how one can compute reduced forms and reduction matrices using 
the Kovacic algorithm for Liouvillian solutions (and its extension to algebraic 
solutions in ISinger and Ulmer although it will probably not surprise 



specialists, we include it for completeness and because the reduction matrices 
will be used in the next section. 

Proposition 7. Consider a 2 x 2 linear differential system Y' = NY with 
Tr(N) = (i.e N E sl(2,k)) and Galois group G C SL(2,C). Then G C 
SL(2,C) if and only if there exists an algebraic extension k° of k such that one 
of the following (mutually exclusive and to be read in this order) cases holds. 
The fi eld k° is given by the Kov acic algorithm and its extension by Singer and 
Ulmer lSinaer and Ulmei\ X199& ) (Theorem 4.1) for algebraic solutions. 

Case (1): The system [N] admits two solution^ Yi,Y% E {k°) 2 . Then Gm is finite 
and q = {0}. 

Let P = (Yi,!^) G SL(2,k°) (after multiplying Yi by a scalar so that 
dct(P) = I). 

Then P[N] =06 fl(fc°) o,nd P is a reduction matrix. 

Case (2): The system [N] admits one solution Y\ E (k ) 2 . Then G N = G a - 

Let P = (Yi,F 2 ) E SL(2,k°) (where F 2 E (fc°) 2 is any vector such that 
det(P) = I). 

Then P[N] = [ " o ] G 8a(k°) and P is a reduction matrix. 

Case (3): There exist a E k° and F\,F2 G (k ) 2 such that, if f is a solution of 
f = af. Y\ := f.F\ and Y% := jF 2 are solutions of [N] (exponential over 
k°). Then G° N — G m . 

Let P = {F\,F2) E SL(2,k°) (after multiplying Fi by a scalar so that 
det(P) = 1). Then P[N] = [ ° _ o ] G Q m (k°) and P is a reduction 
matrix. 

Case (4): There exist a E k° and F\ E {k°) 2 such that, if f is a solution of f = 
af, Y\ := f.F± is a solution of [N] (exponential over k°). Then g = 
{[ o -c ] \c,d E C). Let P = {F 1 ,F 2 ) E SL(2,k°) (where F 2 E (k°) 2 is 
any vector such that det(P) = 1). Then P is a reduction matrix. 

2 This phrase is an abuse of language which, precisely means: the space of the solutions of 
[N] that belong to (fc°) 2 is of dimension two. 
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Proof. We review each case (in the proposed order) and study the Galois group 
in each case. Classifications of the algebraic subgroups of SL(2, C) will show 
that these four cases cover all proper algebraic subgroups of SL(2, C). 
Case (1) is the case when the group is finite. Then the proposed P is a fun- 
damental solution matrix. Our hypothesis on N (i.e Tr(N) = 0) implies that 
det(P) is constant. Multiplying one column by a constant, we may secure that 
det(P) = 1. If now U is a fundamental solution of [N] and we write U = PC, 
then C = ; as C is a fundamental matrix for [P[JV]], this confirms that 
P[N}=0. 

In Case (2), G° has an invariant vector Y\ so G° C G a and, since G is not finite, 
(otherwise case (1)), we have G° = G„. Pick a fundamental matrix U with 
determinant 1 and write it as U = PV. As Yy is a solution, the first column of 
V can be chosen to be (1, 0) T . Now, as det(V) = 1, we obtain V = [ J ° ] for 
some v&K.AsV' = P[N].V, the form of P[N] follows. 

In Case (3), we have g = g m . A fundamental solution matrix U can be written 
U = PV with V = [ o f °i ] and the result follows. 

In Case (4), the galois group is G = { [ c c -i ] \c & C* ,d £ C] . Calculations 
are the same as in (2) and (3) and are left to the reader. □ 

Remark 1. 1. In case (1), the field k ° can h ave arbitrarily large degree (in 



the notations of \Sinaer and Ulmei \l99&) . if G = D^ L2 then k° is of 
degree An over k). 

In case (2), G = I ^ Ji k \ k £ {0, .., n — 1}, d £ Cj where fi is an n- 

th root of unity. Then Y x = f.F 1 with Fi £ k 2 and /" £ k. So k° := k(f ) 
is a cyclic extension of degree nofk. 

In case (3), we have G = G m (then k° = k) or G = D^ 2 (then k° is 
quadratic over k). In case (4), k° = k. 

2. Cases (1), (2) and (3) are the cases when g is abelian. In that case, we just 
saw that the field k° may have large degree. However, for the applications 
that we have in mind (normal variational of Hamiltonian systems), it 
turns out that k° — k in most example^ that we know. 

For 2x2 systems, the construction shows that reduced forms are the most 
natural and the simplest forms into which a system N £ M2{k°) (with qm 
abelian) can be put. The fields over which such r eductions can be perf ormed 
effectively via the Kovacic algorithm are given in lUlmer and Weill (Il996h see 
subsection 4.1). 



3 Actually, in all examples from mechanics that we can think of, but we may have missed 
a few. 
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4 Reduced form of first variational equations 
and non-integrability 



Consider a linear differential system [A] : Y' = AY with A £ sp(4, k). Call 
G its Galois group and g its Lie algebra. We further assume that we know a 
rational solution of [A] (i.e. with coefficients in k) so that the material applies 
to variational equations of Hamiltonian systems along a known non-constant 
solution. 

The aim of this section is to establish an algorithm which either finds an 
obstruction to the abelianity of q or returns a reduced form for A if g is abelian. 
In view of its application to Hamiltonian systems, this would either prove the 
non-integrability of the system (an effective version of Theorem [T]) or put the 
first variational equation in a reduced form so as to simplify the study of higher 
order variational equations. 



4.1 Admissible base fields 

In view of practical computation, we first assume that k is an effective field (i.e 
one can perform the arithmetic operations +,—,*,/ and algorithmically test 
when two elements of k are equal). In order to check effectively the abelianity 
of g, we need three algorithmic tools: 

1. The Koyacic a l gorithm (for solv i ng sec ond order differential systems, see 
Kovacid (|l986l ); iuimer and Weill (| 19961) ). 



2. An algorithm to solve Risch equations : for fig £ k, decide whether the 
equation y' — fy + g has a solution belonging to k (and, if yes, compute 
it), 

3. An algorithm to solve Limited integration problems : for f,g £ k, decide 
whether there is a constant (3 and h £ k such that / + fig = h! (and, if 
yes, compute them). 

We now study conditions on k so that all computations required by the above 
problems can be performed effectively. 

Definition 8. We say that k is an admissible field if it is an effective field and 
if it comes equipped with algorithms to perform the three tasks above: Kovacic 
algorithm, solving Risch equations and solving limited integration problems. 

Given an operator L £ k[d], we say that a solution y of L(y) = is rational 
if y £ k, and we say that y is exponential if y'/y £ k. 

Lemma 9. Let (k,d) be an effective differential field such that, given L £ k[d], 
one can effectively 

1. compute a basis of rational solutions of L{y) =0 (i.e solutions ink) when 
such a basis exists, and 
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2. compute all right-hand first order factors (i.e solutions y such thaty'/y £ 
k, exponential solutions) of second order linear differential equations. 



Then k is an admissible field. 



Proof. If y satisfies a limited integration problem, then y is a solution in k of 
L(y) = with L = LCLM(d - f'/f,d - g'/g)d (where LCLM deno tes the 
least common left multiple in k[d], see Ivan der Put and Singer ( 2003 ) p. 50). 
Indeed, a basis of solutions of L is 1, f /, J g. If now y solves a Risch equation, 
then a similar calculation sho ws that y is a solution in k of L(y) — with 
L = (gd — g'){d — f). Last, in lUlmer and Weill (Il996l) . it is shown that, if the 
two conditions of the Lemma are fulfilled, then one can apply a full Kovacic 
algorithm. □ 



In lSinger (1991), Lemma 3.5, it is shown that if k is an elementary extension 
of C(x) or if k is an algebraic extension of a purely transcendental liouvillian 
extension of C{x), then k satisfies the conditions of the lemma and hence is an 
admissible field. 



Remark 2. In Singer 1 199 ft ). Theorem 4-1, it is shown that an algebraic ex- 
tension of an admissible field is also an admissible field. This is used in the 
next sections, where an algebraic extension of k is (sometimes) needed in order 
to perform reduction. 



4.2 Reduction of the normal variational equation 



Let K denote a Picard-Vessiot extension for [A]. In view of an application 
to variational equations of Hamiltonian systems, we assume that we know a 
rational solution of [A] . This allows us to apply the gauge transformation given 
in (JSJ , yielding an equivalent system of the form 



i-N 








a 13,n 


Ol4,n 





n n 


Gl4,n 


"12 

















"21 


— ai2,n 


-"11 



and N := 



nn 

"21 



"12 

-"ii 



e*t(2,fc) 



Proposition 10. Using the notations adopted in section \3.4\ suppose that Qjy 
is abelian. Then there exists a symplectic change of variable P/v £ Sp(i, k°) 
that puts An into the form B := Pjsr[Ajsr] given below: 



Case q n = {0} : B = a^Mi + a u M 2 + ai 3 M 3 

Case Qn — 9a ■ B = a 12 Ali + OuAfa + ai 3 Af 3 + a 2 A,M a with 

Case g N = g m : B = a 12 M x + auM 2 + ai 3 M 3 + a 22 M m 
with a 22 — E I ' E where E ^ k° . 



a 2 i 



k° 
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where M\ := 



M n 













1 







D 

D 

D 

•1 

and M„ 



Mo : = 





' 








1— 1 










1 





















_ 




























1 

































-1 





M 3 := 



10 









Proof. Keeping the notations of the previous section, we write Kn for the 
Picard-Vessiot extension of Y' = NY. Applying the Kovacic algorithm as in 
proposition [71 we obtain an algebraic extension k° of k and a reduction matrix 
P e SL(2, k°) for [N\. Let P = (pi.j) be this reduction matrix. 
We extend this to a gauge transformation on An'- 



Pn := 





Pn 


P21 





Pl2 



P22 



where 



Pn 
P21 



P12 
P22 



reduces N. 



The classification in the proposition now follows from the fact that, as Pn is 
(by construction) symplectic, we have P/v[-<4] G sp(4, D 



Remark 3. These matrices Mi satisfy simple relations: 



(7) 



MiM 2 = -M 2 Mx = M 3 and MjM,- = otherwise. 

In particular, M? = so exp(JWj) = Id+ Mi. 

In section[2j we have seen that the differential Galois group Gn of the normal 
(variational) equation N is a quotient of the differential Galois group G of the 
whole system. If g is abelian then qn is abelian as well. The converse may not 
hold because a non-abelian group may have abelian quotients. We will now give 
criteria to detect obstructions to abelianity by means of the coefficients of A 
and a reduced form R of A when G° is abelian. 



4.3 Maximal abelian subalgebras and effective abelianity 
conditions 

At this point, if g n is abelian, we see that A £ t)(k°) and g C f), where A and 
f) can be chosen in table Q] below. To find out whether g is abelian, we will 
not need to compute g completely: we only need to show there is an abelian 
maximal subalgebra m C i) such that fj c m. 

Lemma 11. Let g be the Lie algebra of Y' = AY and let f) be a Lie algebra 
such that g C t) and A G i){k°). Then, g is abelian if and only if there is some 
P 6 GL n {k°) such that P[A] G m(fc°) ; where m is a maximal abelian subalgebra 
oft)- 

4 Or by the computation of A := Pjv[A/v] = PZ 1 {-^nPn ~ P'n)> usm g the fact that Pjy is 
symplectic and therefore P^ 1 = — J 1 Pn J is easily computed 
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Table 1: Systems after reduction of the Normal Variational Equation 



Qn 


{0} 


Qa 


Q?n 


A 


012M1 + ai 4 M 2 
+a 13 M 3 


a\ 2 M\ + ai 4 M 2 
+CL13M3 + a 24 M a 


a\ 2 M\ + CL14M2 
+ai 3 M 3 + a 22 M m 




spanc(M u M 2 ,M 3 ) 


span c (M a ,M 1: M 2 ,M 3 ) 


span c (M m ,Mi,M 2 ,M 3 ) 



Proof. Suppose q is abelian. Then there will exist some maximal abelian sub- 
algebra m such that q C m. By Theorem [3l we know that there exists P G 
GL n (k°) such that P[A] G C m(fc°), and we are done. The converse 

follows from Theorem [3] and its corollary. □ 

Lemma 12. Keeping the notations from tableUl we have: 

1- if Qn — {0}, the maximal abelian subalgebras oft) are of the form span c (aiMi + 
a 2 M 2 , M 3 ) where (ai , a 2 ) G C 2 . 

■2- */flAf = Qa> fie maximal subalgebras are span c (M 2 , M 3 , M a ) and span c (M a + 
a x Mi + a 2 M 2 , M 3 ) wift («i , a 2 ) € C 2 . 

-tffliv = Qm, the maximal abelian subalgebras oft) are oftheformspan c (M m + 
aiMi + a 2 M 2 , M 3 ) with (ai , a 2 ) G C 2 . 



Proof. The results in the Lemma are easily deduced from the multiplication 
tables given below 



[Case Q N = Qa] 


[Case W = Q m } 




M 


Ma 


Mi 


M 2 


M 3 






[,} 


M m 


Mi 


M 2 


M 3 




M a 





-M 2 








M m 





-Mi 


M 2 





Mi 


M 2 





2M 3 





Mi 


Mi 





2M 3 





M 2 





-2M 3 








M 2 


-M a 


-2M 3 








M 3 














M 3 















1. If = 0, in the table we see that [Mi , M 2 ] = 2M 3 7^ whereas 
[Mi, M 3 ] = for i G {1,2}. Therefore, any abelian subalgebra of \) 
must necessarily be a subalgebra of an algebra of the type span c (aiAfi + 
a 2 M 2 , M 3 ) where (a% , a 2 ) G C. 

2. If 0at = Q a , in the table we see that [Mi , M 2 ] = 2M 3 ^ and [Mi , M a ] = 
-M 2 j£ whereas [Mi , M 3 ] = for i G {1,2}. Therefore, any abelian 
subalgebra of [} containing M a must necessarily be a subalgebra of 
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- either an algebra of the type span c (M a + a\M\ + a 2 M 2 , M3) where 
(ai , a 2 ) € C. 

- or span c (M a , M 2 , M 3 ) where (ax , a 2 ) G C. 

3. If qn = Qm- since the algebra span c (M m , Mi , M 2 ) is not abelian and 
has no abelian subalgebras (other than the monogenous ones), all abelian 
subalgebras of f) containing M m must be of the type span c (M m + a\M\ + 
a 2 M 2 , M 3 ) where (ai , a 2 ) € C 2 . 

□ 

Remark 4. In the above lemma, two cases may be simplified. In the second 
part of case (2), span c (M a + ot\M\ + a 2 M 2 , M3) is conjugate to span c (M a + 
aiMi , M 3 ) (via P = lA + a 2 M x ). In case (3), span c (M m + aiMi+a 2 M 2 , M 3 ) 
is conjugate to span c (M m , M3) ('nia P = Id + ct\M\ — a 2 M 2 ). 

Before we proceed to our main result, we need a lemma about the structure 
of solution matrices. 

Lemma 13. Let Y' = AY with A € sp(4,fc°) as given in table [71 vaae \15[ 
Assume that is abelian. Then, depending on Qn , a fundamental (symplectic) 
solution matrix U £ Sp(A, K) can be chosen in the table below: 



0N 


aid = {0} 


Qa 


0m 


u 


r 1 Qi q 3 n 2 1 

1 ^2 

10 
[_ -fii 1 J 


I" 1 fii fl a Q 2 + LQ! 1 

1 n 2 l 

001 

[ -Qi 1 J 

L ■- fa 2 4 i k° 
with a 2 4 £ k° 


[ 1 EQ,i fi 3 % ] 

e n 2 

10 
L -| J 

£ := exp(J* a 22 ) ^ fc° 
wzt/i a 2 2 G fc° 



Proof. This table is a direct consequence of Proposition[T5]given in the appendix 
(pageHU). □ 



We now state our main result. This next theorem gives a reduction algorithm 
which puts, if possible, A into reduced form. In what follows, the term "P is a 
reduction matrix" means that Lie(P[A]) is abelian (so, strictly, the system may 
be only partially reduced but this is enough for us). 

Theorem 14. Consider a differential system Y' = AY , where A is chosen in 
table [JJ page 1 1 51 Let q be the Lie algebra of its differential Galois group. Then, 
Q is abelian if and only if 

1. Case gjv = {0}: one of the two assertions below holds. 

1. There exist yi,y 2 € k° such that y[ = an and y' 2 = 014. In that case 
q C span c (Ma) and P := Id + y\M\ + y 2 M 2 is a reduction matrix. 
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2. There exist (a\ , a%) G C 2 \ (0, 0) such that the equation y[ = a 2 a\ 2 — 
a\ai4 has a solution y% € k° . In that case g C span c (a2-Mi + 
ol\M-2. , A/3) and P := Id + 2ct * Q2 2/i (aiMj — a^A^) is a reduction 
matrix. 

2. Case g^r = g Q : one of the two assertions below holds. 

1. There exists y\ G k° such thaty[ — avz- In that caseg C span c (A/2, M 3 , M a ) 
and P := Id + y\Mi is a reduction matrix. 

2. There exist (ai , a 2 ) G C 2 \ (0, 0) and 2/1,2/2 G k° such that 

y'l = ai2 - aia 24 ,„\ 
y' 2 = a i4 - a 24 2/i - a 2 a 2 4- 

In that case, g C span c (M a + ct\M\ + a 2 M 2 , M 3 ) and P := Id + 
y\M\ + y 2 M 2 is a reduction matrix. 

3. Case g^v = 0m : system of Risch equations 

V'l = -^222/1 + Ctl2 /g\ 
2/2 = «222/2 + ai4 

/ias a non-trivial solution in (k°) 2 . In that case g C span c (M m , M 3 ) and 
P := Id + 2/1 Afi + 2/2 A-^2 *s a reduction matrix. 

Proof. (Theorem) We first show the case qn = {0} in detail, together with a 
reduction procedure. The next cases will then be easier to follow. 

1. Case qn = 0: In this case (by lemma [T3|) . we have A — a\ 2 M\ + a\^M 2 + 
013 A/3. Lemma [TBI shows that the system Y' = AY has a fundamental 
solution matrix of the form 

U = Id + fiiMi + n 2 M 2 + n 3 M 3 

Assume now that g is abelian. By Lemma [TT] This is equivalent to 

G m Ql:Q2 := span c (aiA/i + a2Af2 , M3) with (ai , 012) G C 2 . 

We pick a G exp(g) C G° and study its action on U. Using relations ((TJ) 
and abelianity of m QljQ2 , we see that a = Id + (3i{aiMi + a 2 M 2 ) + (3 3 M 3 
for some (3i,/3 3 G C. 

By a quick calculation using relations ([7])) the action of a over U is given 
by a(U) = U.a = 

Id + (fii + /3iai)A/i + (il 2 + /3 1 a 2 )M 2 + (fi 3 + f3 3 + /3i(a 2 fii - ai n 2 ))M 3 . 
Identifying this with a applied to coefficients of U, we infer that 

o-(n x ) = n 1 + a 1 /3 1 

a(n 2 ) = fl 2 + a 2 l3 1 



17 



If ct\ = o- 2 =0: this shows that f^i and fi 2 are fixed by exp(g) hence by 



G° (e.g iFulton and Harris! (|1991l) . proof of proposition 8.33) so they are 
algebraic. Now we also have ft^ £ k° hence £ k° (for i = 1, 2). 
We introduce the reduction matrix P := Id + Q^Mi + fi 2 M 2 £ GL^k ) 
and easily check that U = P. (Id + Hence, the new system has 

matrix P[A] = (a 13 + a 12 fl2 — auVL^Mj, and g C span c (M 3 ). 

If now (0:1,0:2) 7^ (0,0), we let / := o 2 f2i — oif2 2 . It is left fixed by G° 
and /' = o 2 ai 2 — 01014 £ k° so / £ k° as above. We use the reduction 
matrix P := Id + 2^7/ («i^i - a 2 M 2 ) € GL 4 (k°) and obtain 

P[A] =^(aiMi+O2M2) + (ai3"2//3)M 3 Gm Ql , Q2 (fc ) where /3 = 
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so g C Tn QljCt2 and the system is again (partially) reduced with an abelian 
Lie algebra. 

Conversely, if one of the conditions of the theorem holds, the reduction 
(gauge transformation) matrices P given above send A to a fc°-point of 
span c (Af3) or m ai ^ a2 respectively, showing that g is abelian. 

2. Cases Qn = 9a- In this case g must always contain M a otherwise we would 
be in the case g n = 0. By lemma [T3l we have A = a^Mi + C114M2 + 
(I13M3 + a2iM a . The above lemma shows that the system Y' = AY has a 
fundamental matrix of the form 

U = Id + fiiMi + ^ 2 M 2 + ^M 3 + LM a - L0iMiM o . 

Assume now that g is abelian. This is equivalent (Lemma fTTT) to either 
C m := span c (M 2 , M 3 , M a ) or g c m QljQ2 := span c (M a + oiM 2 + 
a 2 Mi, M 3 ): 



(a) g C m := span c (M 2 , M 3 , M a ): pick a £ exp(g) C G° and study 
its action on U. Using relations |7]) and abelianity of m, we see that 
cr = Id + /3 2 M 2 + /3 3 M 3 + /3 M a for some /3 2 , #3, /3 £ C. Computing 
the action of a over £/ we obtain Ua = 

ld+n 1 M 1 + (p2+Pin 1 +n 3 )M 3 +{n 2 +f32)^h+Wn 1 -Ln 1 )M 1 M a +(f3+L)M a 

Identifying this with a applied to coefficients of [/, we infer that 
a(fli) = fix which is equivalent to Jai 2 =: fii £ k°. The reduction 
matrix obtained this way is 

P = Id + O1M1 

and P[A] = (ai 3 — 2au fl± + a 2 4 Clf)Mi + (au - a2A ^i)Af 2 +a 2 4 M a . 
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(b) fl C m := span c (aiAfi + a 2 M 2 + M a , M 3 ): pick a £ exp(g) C G° 
with a = exp(/3i(M a + a x Mx + a 2 M 2 ) + /3 3 M 3 ). We compute U.a 
and compare this with entries of a(U) to obtain: 

(ll + anpi 
L + fa 

= n 2 + q 2 /3i - Lp x ai - %i 




Let J/i := fii — Qfii; the first two relations show that y\ is fixed by 
all <7. Hence yi € k° and j/i is a solution in k° of t/J = a\ 2 — a\a 2 ^. 
Similarly (but with a few lines of calculation), letting y 2 :— il 2 — 
a 2 L + ^a\L 2 (i.e y' 2 — 014 — a 2 ^y\ — a 2 a 2 4), the third condition 
shows that y 2 £ k° The reduction matrix obtained this way is 



P m = ld+yiMi+y 2 M 2 and 
with h = y' 2 y x - y[y 2 + a 13 . 



P m [A] = a2i(M a +aiM 2 +a 2 Mi)+hM 3 



3. Case g^r = Q m '■ Again, M m must belong to q (otherwise we would be in 
the case Qn = 0). By Lemma [TBI we have A = a 12 M x + a 14 M 2 + a 13 M 3 + 
a 22 M m . Lemma [TBI shows that the system Y' = AY has a fundamental 
matrix of the form 



U := 



" 1 


EQx 


fi 3 


Q.2 ' 

E 





E 


n 2 











1 













1 

E . 



Assume now that q is abelian: by Lemma this is equivalent to g C 
span c (aiMi + a 2 M 2 + Af m ,M 3 ) for some ai,a 2 £ C. We pick a a £ 
exp(m) C G°: notice that since 7^ the form of a will be 



M„ 



1 ai (e c — 1) oi\ a 2 (e c — e c — 2c) + (3 a 2 (1 — e 
e c a 2 (e c -l) 
1 
-ai(e c -l) e" c 



5 ) 



with c £ C*,/3,ct!j € C for i = 1,2. Computing cr(J7) = U.M a we obtain 
the following relations linking E, Qi, Q 2 and the constants linked to a: 



a(E) 
<r(fii) 

£7(fi 2 ) 



e c .E 

«l 1 ai 

Sl 2 + a 2 c c E — Ea 2 . 



Again, it is easily checked that £7i — ^ and f2 2 — a 2 -E are fixed by all a. 
This is equivalent with saying that the following equations 



Vi = 
V2 = 



-a 22 yi + <zi2 
a 22 y 2 + an 
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have solutions 2/1,2/2 € k°. The reduction matrix obtained this way is 
P m = ld+y 1 M 1 +y 2 M 2 and P m [A] = (a 12 y2+a 13 +a 14 y 1 )M 3 +a 2 2M m . 

□ 

Remark 5. In the case qn — Qm, we actually proved that q is a subalgebra of 
span c (M3 , M m ) instead of the maximal abelian subalgebra m = span c (aiMi + 
a 2 M 2 + M m , M3). Indeed, if we pick a differential system [B] Y' — BY with 
B := f(x) («iMi + a 2 M 2 + M m ) + g(x)M 3 G m(k) the abelianity conditions 
given in the theorem are written as: 

f y'i = -/OX yi - an ) 
1 y'2 = f( x )( V2 + a 2 ) 

and this system always has a non-trivial (constant) rational solution (jji , y 2 ) — 
(cti , —a 2 ). Therefore, there exists a linear change of variables given by P m := 
Id + ol\M\ — a 2 M 2 such that P m [B] — g(x)Ms + f(x)M m . This recovers what 
we saw in remark^ the fact that these two algebras are conjugate. A similar 
thing (using again remark^ holds for the second additive case. 

4.4 A reduction algorithm 

We start with a partial reduction algorithm which summarizes the procedures 
established above. 

INPUT: A matrix A 6 sp(4, k), with k an admissible differential field, and a 
particular solution of Y 1 = AY. 

Step 1: Perform symplectic reduction using section 2. 

Step 2: Reduce the NVE applying Kovacic's algorithm as in section 2 
(gives the reduction field k°). 

Step 3: If Lie(NVE) is abelian, follow the algorithm of Theorem [TJ] 

OUTPUT: either "g non-abelian" or P G GL(4, k°) such that Lie(P[A]) is 
abelian 

We call B := P[A] the output of this algorithm and let m := Lie(B). We 
could now wish to complete the reduction in the case when q C m. The Lie 
algebra q is a subalgebra of m but not any algebra. Indeed, when qn is Q a 
(resp. m ), the matrix M a (resp. M m ) must belong to q. Otherwise we would 
fall again in the case 0at = 0. We give in the next lemma the list of possible g. 

Lemma 15. List of all possible abelian subalgebras: 

1. The subalgebras of m — span c (aiMi + a 2 M 2 , A/3) are among {0}, m, of 
the form span c (aiAfi + a 2 M 2 ) with (ai , a 2 ) G C 2 \(0,0), or span c (M 3 ). 
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2. The subalgebras of m = span c (M + ct\M\ + a 2 M 2 , M3) that contain M a 
are among {0}, m, or of the form 8pan c (aiMi + a^M?, + (X3M3 + M a ) 
with (a 1 , a 2 , a 3 ) G C 3 . 

5. T/ie subalgebras ofm — span c (M2 , -M3 , M a ) that contain M a are among 
{0}, m, of the form span c (a 2 M 2 + a 3 M 3 , M a ) with (a 2 , 0^3) G C 2 , 
or 0/ the form span c (M 2 , a 3 M 3 + M a ) with a 3 G C* . or of the form 
span c (M 3 , a 2 M 2 + M a ) with a 2 G C*, or span c (a 2 A/ 2 + a 3 M 3 + M a ) 
with (a 2 , 03) G C 2 . 

4. The subalgebras ofm — span c (M m , M3) that contain M m are among {0} ; 
m, or of the form span c (Af m + a 3 M 3 ) with a 3 G C. 

Proof. We leave the computation to the reader. □ 

We may infer a reduction algorithm from the a bove lemma. Name l y, the 
proofs of propositions 5, 6, 7 and 8 in chapter 4 of lAparicio-Monfortei ( 2010l ) 



(one proposition per numeric item of Lemma ll5j) sum up a reduction algorithm. 
Furthermore, the reduction will be maximal because the procedure proposed 
deals only with abelian Lie algebras. An outline of the proof of those four 
proposition runs as follows. We pick for each case an element a G exp(g) 
and compute its action over a symplectic fundamental matrix of solutions U. 
From the relations obtained by means of these computations we extract the 
reducibility conditions. The reduction matrix and corresponding reduced form 
are obtained using a direct computation. To prove that P [B] is reduced instead 
of just partially reduced, we prove either that g is monogenous or that dimc(g) 
is equal to the minimal number of generators of g as a Lie algebra. At the 
outcome, we obtain a fully reduced form for our differential system. 

Remark 6. When a system is in reduced form, one may sometimes yet find 
a simpler reduced form. For example, assume we were able to write some co- 
efficients as a = f + g with f,g G k° ; for instance, if k = C(x) we can use 
partial fraction decomposition and Hermite reduction to obtain a = f + g ( and 
g with simple poles). Then, like in the reduction process, we can find a gauge 
transformation which will remove the f part of a, leaving only g. 

Example 16. Consider the systemY' = MY given by M = ]d ^-M 1 + x3 ~^l +1 M 3 G 

Sp(2,C(x)) with f\(x) — and f 3 (x) — x3 ~^l +1 ■ Its Lie algebra is g M = 

span c (M! , M 3 ). It is easily checked that M is already in a reduced form. Using 
partial fraction decomposition, we obtain: fi(x) = and f 3 {x) = x 2 + — 

Now, applying the gauge transformation Q := Id + \M\ + \M 3 G GL4(C(i)) 
we obtain the equivalent and simpler system Q\M\ — —M\ + -^M 3 . 

This kind of procedure can be generalized to any monomial and e lementary 



extens ions of k where k is a differential held of characteristic 0, see iBronstein 
( 20051) . 
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5 Reduction and integrability for Hamiltonian 
systems in dimension n = 4 



5.1 The algorithm 

The reduction algorithm shown in the previous section is easily turned into an 
effective version of the Morales-Ramis criterion for meromorphic Hamiltonian 
systems with two degrees of freedom. Indeed, consider a Hamiltonian system 
such as (JJ) with two degrees of freedom and assume that Y is one of its non- 
constant integral curves. The following algorithm either returns a reduced form 
of the variational equation of (JXJ) along Y, or tells us that the system JT]) is not 
meromorphically Liouville integrable. 

INPUT : A meromorphic Hamiltonian function H : U C C 2 ™ — > C, a partic- 
ular non-stationary integral curve z(t) for the system ([1]). 

1. Apply the abelianity test to A :— Hess(H)(z(t)). 

2. If fj is abelian apply the reduction algorithm otherwise stop. 

OUTPUT : If g is abelian we obtain a reduced form R £ g(k°). Otherwise g 
is not abelian and the Hamiltonian system (JT]) is not meromorphically Liouville 
integrable. 



5.2 An application: the lunar Hill Problem 

Let us now apply our algorithm to a real problem: a simplification of the Sun- 
Earth-Moon problem, more widely known as the Hill Problem. Meromorphic 
n on-integrability of this syst em has been proved by Morales, Ramis and Simo 
in lMorales-Ruiz et al. I (l2005h . Applying our techniques, we give a new (simpler) 



proof of this result. Although, the lunar Hill problem is modeled by a Hamil- 
tonian functio n with three degrees of fr eedom like all three bodies problems, 



it is shown in iMorales-Ruiz et al.1 (j2005l ) that via a symplectic change of vari- 



able over fc, the lunar Hill Problem can be seen as the two degrees of freedom 
Hamiltonian system given by the Hamiltonian function: 

H := i(qiq 2 - p\p 2 ) - ^q\q 2 (q\P\ - 92^2) - 4i(3gf - 2qfqj + iql)qiq 2 . 



We compute Xh (see Morales-Ruiz et al. ( 20051 )) and notice that it possesses 
two invariant planes: 

Hi := {q 2 = , pi = 0} and Yl 2 := { qi = , p 2 = 0} . 

One checks easily that Xh has a particular solution over the invariant plane Hi 
given by 7(35) = (/(at), 0, 0, if'{x)) such that f{x) satisfies (f'{x)) 2 = -f(x) 2 + 
4/(x) 6 + 2h so that, by differentiating, f"{x) = -f(x) + I2f(x) 5 . We may 
take f(x) 2 — 3 p(v )+1 where h G C is an indeterminate constant and p(x) = 
p(x; |, Jp(l — 216/i 2 )) is a Weierstrass-p function. We work over the differential 
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field k := C(f(x),f'(x)). We refer to the elements of k as rational functions. 
Consider now the variational system over 111 along 7. Its matrix A is 



.4 



-4/(x) 2 

-i 

-i(l - 60/(x) 4 ) 

-i(l - 60/(:r) 4 ) -8i/(as)/'(x) 4f(x) 2 



We perform step 1. of the algorithm: we compute P and An ■= P[A] 
P- 1 {AP-P I ): 





\ /'(*) 








" 







1 








P := 





/'(*) 


1 

/'(*) 















1 . 



.4 



" /'(*) 

/'(*) 



F(z) 







/'(*) 


4/(^) 2 



/'(*) 






where F(x) := hpj^ '- 

N is defined by: 



N 



Thus, the normal variational equation matrix 



F(x) 







Variation of constants shows that the Normal Variational Equation has a 
fundamental solution matrix 



U N := 



/'(*) 




1 

77?) 



Thus we have Un G Sp(2, k) which is equivalent (in the notations of proposition 
[7] on page fTUj) to qm — 0. Now reduce the normal variational equation: the 
matrix Pjv G >?P(4, fc) puts j4at into the form A := Pn[An] where 



Pn := 





/'(*) 


8i{f(xf-hf(x) 2 






1 

- 









and A := 



4G(s) 












-4G(x) 



with G(z) = 4 f{x) ^ f{x) 2 +2h — L G fc. 

We apply our reduction algorithm (theorem [TH page ITS)) . We are in case 1 
so we search for (non-zero) constants (a, j3) £ C 2 \ (0 , 0) such that J (a4G(x) + 
P e ^- Changing variable xtot — f(x)mI (so = -^=5=5==), we are 

reduced to a limited integration problem on a hyperelliptic curve. Using section 
4.1 and setting h = 1, this is in turn equivalent with finding (a , 0) E C 2 \ (0, 0) 
such that 

I := J a4Fi(t) + (3F 2 (t)dt e C(i, \J \& -t 2 + 2) (10) 



7W 
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where 



4t 2 (-2t 6 + t 2 -4) 



(4i 6 - i 2 + 2) 3 / 2 ' w ' (4i 6 - t 2 + 2) 3 / 2 ' 

To solve this algorithmically, we apply our method from lemma l9l page [T2l Each 
Fi is a solution of Li :— ^ — ^|jy for j = 1,2. Hence, ai<i(i) + (3F 2 (t) will be a 

solution of LCLM(L\ , L2) for all (a, /3) £ C 2 which implies that I is a solution 
of the differential operator L := LCLM{L\ , £2)37- Thus, solving our limited 
integration problem will amount to finding non-trivial rational solutions (i.e 
solutions in C(£, V 4i 6 — t 2 + 2) other than 1) for the (parameter free) differential 
operator 

_(d\ 3 (44i 6 -3i 2 -2) /d\ 2 t 2 (48t 8 -24i 4 + 96t 2 - l) d 
L: _ W + (-t 2 +4t6 + 2)t (.dij +3 (-t 2 + U^ + 2) 2 

By studying the singularities of L, we prove that the latter does not have non- 
trivial rational solutions. Indeed, the equation has three regular singularities: 
with exponents {0,1,3}, a such that 4a 6 — a 2 + 2 = with exponents 
{0 , 1/2, 3/2} and last 00 with exponents {0,0, 8}. All these singularities 
have minimal exp onent equal to s o any solution of L in C(i, -\/4i 6 — t 2 + 2) 
must be constant (Sin ger and Ulmer 

Therefore, Theorem [TJ] shows that g is not abelian which implies that the 
Hill Problem is not meromorphically integrable. 

Remark 7. For the meromorphical non-integrability of the Hill problem, one of 
the referees suggested an alternate proof using integration on an algebraic curve 
instead of the language of differential operators. Let us give its outline. Taking 
as above h — 1, we can write the integral I as 

I := / P ® dt 



D(t)y/m$ 

where P(t) := 4at 2 (-2x 6 + x 2 - 4) + i(3 and D(t) := At 6 - t 2 + 2. Proving our 
point, amounts to proving that unless we pick (a , f3) = (0 , 0), the integral I ^ 
C(t, y/D(i)). We assume that if I is rational then there exist some A(t), B(t) € 
C(t) such that I = A(t) + B(t)y/ D(t) . Differentiating this expression, we obtain 
A'(t) — and therefore 

B >(t) + m B{t) - 41 = 0. (id 

1 ' 2D(t) y ' D(t) 2 K 1 

Solving this Risch equation we are done. Indeed, we see that the only rational 
solution it admits corresponds to [a. , f3) — (0, 0) and is trivial. It is true that 
this proof only requires the resolution of a Risch equation instead of handling 
a third order differential operator. However, we think that our first choice is a 
good one. Indeed, we avoid discussing parameters (a, /3) and our method can be 
applied systematically to any situation satisfying the conditions stated in section 

4.1. 
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Remark 8. 1. if we had chosen the energy level h = 0, the same argument 
leads to looking for solutions in C(t, V At e — i 2 ) of the differential operator 

(d_\ 3 (44^-3) / d\ 2 (48f 8 -24f 4 -l) d 

\dt) + t (2*2-1) (2*2 + 1) {&) + t 2 (2t2-l) 2 (2i2 + i) 2 rfi 

and, indeed, we find (using MAPLE ), the solution ^t^u/sit^ it' 1 ^ ^' 
So, on the energy level h = 0, our reduction method shows that the system 
has an abelian Lie algebra. 



2. In an example like this where coefficients are parametrized by Weierstrass 
functio ns, one would need in general to use the spe c ial al gorithms devel- 
oped in \Sinaei\ \l99i ) and improved in \Buraer et all \200A ) to achieve the 



reduction. 



Remark 9. Notice that the normal vari a tiona l equation computed in our proof 
and the one given in \Morales-Ruiz et al. 1 200 A ) are different because they arise 



from two different definitions of Normal Variational equation. On one hand, 
the construction we use can be termed as algebraic normal variational equation 
since it is obtained by a purely algebraic manipulation (see section [2.2.2\) : to wit, 
a linear Hamiltonian change of variable obtained from a parametrization of the 
integral curve T via the symplectic Gram - Schm idt algorithm. On the other hand, 



the notion used in \Morales-Ruiz et al. {200 A ) could be qualified as geometric 



normal variational equation since reduction (in general not symplectic reduction) 
is performed with respect to an invariant plane that contains the nondegenerate 
integral curve T . Both notions coincide in the case when the invariant plane 
considered in the geometric construction is actually a symplectic manifold (with 
symplectic form the restriction to the invariant plane of the global symplectic 
form). 



6 Conclusion 

The notions of reduction and reduced form developed in this paper provide 
a procedure to decide the abelianity of the Lie algebra of the differential Ga- 
lois group of the variational equations [A] of Hamiltonian systems. Previously, 
applications of the Morales-Ramis criterion were generally limited to normal 
variational equations. 



When the Lie algebra is indeed abelian, putting the system into a reduced 
form is very convenient because it also allows to (partially) reduce higher vari- 
ational equations in view of a concrete application of the Morales-Ramis-Simo 
criterion. On one hand, the higher variational equations are reducible linear 
differential systems whose diagonal blocks are A and its symmetric powers (in 
the sense of Lie algebras) ; knowing a (partial) reduction matrix P for [A] , its 
symmetric power Sym m (P) is a (partial) reduction matrix for sym m (A) hence 
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inducing a partial reduction on the higher variational systems. Applying the 
techniques of se ction 4 to develop constru ctive ab elianity criteria for such sys- 
tems is used in Aparicio-Monfortel ( 2010h and in Aparicio-Monforte and Weil 



( 20111 ) and will be the subject of other future work. 



In many cases, Hamiltonian systems come as parametrized systems, for ex- 
ample with a base field C — Q(t%, . . . ,t s ). Though there cannot exist an algo- 
rithm deciding for which values of the p arameters the variational equation will 



admit rational solutions (jBoucherl (|2000l) h it turns out that in many situations, 



authors have been able to apply criteria like the Kovacic algorithm (or variants) 
to overcome that difficulty We note that, in this case, the problem of apply- 
ing our techniques becomes tractable. As shown in section 4, when the normal 
variational equation has an abelian Lie algebra, the abelianity of g depends on 
whether integrals belong to k, the latter depending on whether some residues 
are null or not: this problem should be decidable. So, we believe that for fami- 
lies of parametrized Hamiltonian systems, once the normal variational equation 
has been fully reduced (which is not the contribution of this paper) , the remain- 
ing part of the reduction should be tractable even in the presence of parameters. 



Appendix: Symplectic Gram- Schmidt Algorithm 
and Symplectic Linear Differential Systems 

The material in this appendix is mostly well-known and included for the sake 
of the exposition's clarity. 

A.l A symplectic Gram-Schmidt method 

Let (V,lu) denote a symplectic vector space of dimension 2n with a basis u := 
{ui,. . . ,U2n}- We briefly review a stan dard construction sometimes ca lled sym- 
plectic Gram-Schmidt method (see e.g I Abraham and Marsdenl ( 19781 ). chapter 



3.1) for computing a symplectic (or Darboux) basis e := {ei,. . . ,e2 n } from u, 
i.e one on which the matrix of u> is J. 

Let e\ := U\. As u is non-degenerate, one of the u,, say u 2 satisfies 
w{e\,Ui) 7^ hence we may set e, i+1 := u ^ u 2 so that w(ei,e n +i) = 1. 

Let Vi :— span c {ei, e„+i} and V2 ■= V x " be its symplectic orthogonal. Then 
Vi p| V2 = {0} and a basis of V2 is given by v%-2 '■= Ui — uj(v,i, e n +i)ei + 
uj(m, ei)e n +i (for i — 3,...,2n) so V = Vi ® V2 and we may apply recur- 
sively the above to the basis vi, . . . , v n -2 of V 2 - By construction, the result is a 
symplectic basis. 

A. 2 Symplectic linear differential systems 

Let [A] : Y' = AY with A € M.2n{k) be a linear differential system such that 
Gal([A}) C Sp(2n,C). We may transform A into a gauge-equivalent matrix 
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B 6 sp(2n, C) as follows. 

To [A] is associated a differential module (M. , d) where M — k 2n and the 
action of a basis and the action of d over a basis e = {ei, . . . , e 2 „} is given 
by d(e) = -*Ae. This way, for F = ^V^e* £ M, 8Y = Y' - AY. As 
G C Sp(2n, C), there exists a vector w — (. . . , itfjj, . . .) £ /\ 2 M (the coordi- 
nates lOjj € fc are expressed in the basis (e^ A €j)) such that 9(w) = 0. The 
matrix il — {uJi,j) is a skew-symmetric non degenerate matrix (with the conven- 
tion Wi t i = 0). This n is the matrix of a symplectic form u on M. Applying the 
above symplectic Gram-Schmidt method, a symplectic basis e = {e\, . . . , e„} of 
M. is obtained. Let B be the matrix of d on this basis. As the matrix of u on 
e is J, the matrix of the (musical) isomorphism ur between A4 and A4* is still 
J so that JB + f BJ = 0. 



The construction of Normal Variational Equations in section 12.2.21 admits 
a simple explanation in this formalism. We assume that A £ sp(2n, C) so 
that the matrix of the symplectic form on M is J in this basis. If we know a 
particular solution Y\ = t (yx, ■ ■ ■ , ysn) £ k 2n of Y 1 = AY, we obtain what we 
call the (algebraic) Normal Variational Equation by applying the Gram-Schmidt 
method to the basis Y\ , e 2 , . . . , &2n of Ai . Notice that the Normal Variational 
Equation defined in this alge braic way doesn't necessa rily coincide with the 
definition given for instance in Morales-Ruiz et al. 

In our work, the construction implies that y n +i is by construction a constant 
vector : by duality, J ■ *(1, , - for. 1 , ) is a s olution of V ; = — * AY and therefore 
is a first integral of V = AY (|Weill (|l995[ ); iMorales Ruij (jl999t) ). 



A. 3 Symplectic solution matrices of symplectic linear dif- 
ferential systems 

We consider again [A] : Y' = AY with now A £ sp(2n , C) and recall how to 
construct a symplectic fundam e ntal solution matrix (this result is proved e.g in 
Ivan der Put and Singerl (|2003h : iMitschi and Singer] (|2002h V 



Lemma 17. Let k be a differential field and let [A] : Y' = AY with A £ 
Sp(2n, C) be a linear Hamiltonian system. Let K D k be the Picard-Vessiot 
extension of k associated to [A] . Then, [A] admits a symplectic fundamental 
matrix of solutions U £ Sp(2n, K). 

Proof. Assume that we have a fundamental solution matrix U . Then JU is a 
fundamental matrix of the dual system [A*] : Y' = —*AY so there exists a 
constant matrix $ such that JU = *C/ _1 $. Therefore, $ = t V JU which is 
antisymmetric so it is the matrix of a symplectic form over the solution space of 
[A] spanned by the columns U±, . . . , U2n of U. Applying the symplectic Gram 
Schmidt method to the basis Ui, . . . , Uin (for the symplectic form uj(Ui , Uf) := 
t U^>Uj) yields a symplectic basis Vi, . . . , Vm and the matrix V = (Vi, . . . , Vin) 
is a fundamental solution matrix which satisfies by construction, t VJV — J. □ 
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Remark 10. As noted by a referee, there is a simpler proof when k = C(x) 
or k is some field of germs of holomorphic functions. Assume that x = 
is an ordinary point and construct the local fundamental solution matrix U £ 
A^ n (C[[a;]]) such that U{Q) = Id. Then, as l AJ + J A — 0, a calculation shows 
that (tUJUy = and, as U(0) = Id, l UJU = J. 

Now we focus on the case where 2n — 4 and study the structure of symplectic 
matrices belonging to Sp(4, k) := { U £ Mi(k) : t UJU — J}. The Lie algebra 
sp(4, k) is the set of 4 x 4 matrices A £ M^k) such that 



.4 



M 
Si 



s 2 

-*M 



with M,Si€M 2 {k) and t S i = S i 



Consider the linear differential system [A] : X 1 — AX with A £ Sp(4,k). 
Suppose that it admits at least one rational solution X\ € fc 4 . Pick P € 
GL(2n, k) such that X = PY and Pe x = X x where e x := '(1,0,0,0). This 
implies that Y\ — e\ is a solution of the gauge equivalent system Hence, 
the first column of P[A] is null. If, in addition, A £ 5p(4, fc) and P £ Sp(4, k), 
then P[A] G sp(4, k) and its form is 



P[A] = 



a 
nn 


nil 



d 
e 


—a 



e 

nii 


Tin 



(12) 



Proposition 18. There is a symplectic fundamental matrix of solutions, U £ 
Sp(4, K) of the partially reduced system (j which has the following structure 



u = ( u i,j) 



1 




n :i 


912^1 + 





9n 


n 2 


912 








i 








912 


-fi a 


922 



922^2 



where Q :— (qi.j) is a unimodular fundamental matrix of solutions of the system 
Y' = NY, with N := (n itj ) £ sp(2 , k). 

Proof. We can set without any loss of generality that the first column of U :— 
(ui t j), denoted U x , is e x . Since U £ Sp(4,if) (K being the Picard Vessiot 
extension), we obtain S^j = u>(U±, Uj) = u>(e±, Uj) — u^j. Therefore, the matrix 
U can be written in the form 



U = 



1 x fi 3 y 

qn f2 2 912 

10 

912 -^1 922 



and the expression of x and y in terms of the remaining coefficients follows from 
the relations aj(U~2,U3) = ui(Us,Ui) = 0. Furthermore, since uj(U~2,U4) = 1 = 
det(Q), Q is unimodular and we are done. □ 
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